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Abstract 



We study quantum mechanical systems with "spin" -related contact in- 
teractions in one dimension. The boundary conditions describing the 
contact interactions are dependent on the spin states of the particles. 
In particular we investigate the integrability of A^-body systems with 
5-interactions and point spin couplings. Bethe ansatz solutions, bound 
states and scattering matrices are explicitly given. The cases of gener- 
alized separated boundary condition and some Hamiltonian operators 
corresponding to special spin related boundary conditions are also dis- 
cussed. 



Quantum mechanical solvable models describing a particle moving in a local singu- 
lar potential concentrated at one or a discrete number of points have been extensively 



discussed in the literature, see e.g. P, |g, [Tj] and references therein. One dimensional 



problems with contact interactions at, say, the origin (x = 0) can be characterized by 
the boundary conditions imposed on the (scalar) wave function ip at x = 0. The history 
of this problem is well described in 0, §]. It was suggested to divide these conditions 
into two disjoint families: separated and nonseparated boundary conditions, correspond- 
ing to the cases where the perturbed operator is equal to the orthogonal sum of two 
self-adjoint operators in L2(— oo,0] and L2[0, oo) and when this representation is im- 
possible, respectively. Classification of one dimensional point interactions in terms of 
singular perturbations is given in |jl9|. In the present paper we are interested in model 
few-body problems with pairwise interactions given by such potentials. The first model 
of this type with the pairwise interactions determined by delta functions was suggested 
and investigated by J.B. McGuire and C.A. Hurst |2^, ^ |26|, ^ |28|. The eigenfunctions 
for the system of identical particles interacting via delta potentials are given by Bethe 
Ansatz. Intensive studies of this model applied to statistical mechanics (particles having 
boson or fermion statistics) by C.N.Yang and his collaborators lead to the famous Yang- 
Baxter equation 0, |31|. It has been shown in |20|, [18[ that A^-particle systems 
with three-body interactions do not have eigenfunctions given by Bethe Ansatz. In 
the integrability of one dimensional systems of identical particles with general con- 
tact interactions described by the boundary conditions imposed on the wave function was 
investigated. It was shown that the A^-particle system satisfies a Yang-Baxter relation 
not only in the ^-interaction case, but also for two other one parameter (sub)families, one 
with nonseparated boundary conditions and another with separated boundary conditions. 
This fact is not surprising, since the Yang-Baxter equation has been derived for particles 
with boson or fermion statistics. Suppose that the system of A^ particles satisfies one 
of these statistics. Then the eigenfunction equation can be reduced to an equation in 
the sector xi < X2 < . . . < xn, since the value of the total wave function in the whole 
space can be reconstructed using symmetry properties of this wave function. The 
boundary conditions on the total wave function are transferred into certain conditions at 
the boundaries of the sector for the reduced wave function. In fact all three families of 
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boundary conditions obtained in @ correspond to one reduced problem. Hence as far 
as particles with statistics are concerned the only difference between the three families 
is due to the symmetry properties of the wave function, i.e. the rule how the total wave 
function can be reconstructed from the reduced one. Considering particles without any 
statistics the eigenfunctions corresponding to the boundary conditions from the three one 
parameter families the eigenfunctions can be calculated using Bethe Ansatz. In ^ it is 
shown that not only the models satisfying the Yang-Baxter equation have eigenfunctions 
of the type of those constructed following Bethe Ansatz. This is possible, since to derive 
Yang-Baxter equation from Bethe Ansatz one has to use symmetry properties of the wave 
function determined by the statistics. The family of such model operators is described by 
two real parameters. One of these parameters is redundant in the sense that the operators 
corresponding to different values of this parameter are unitary equivalent. It is shown in 
1^ that the redundant parameter can be interpreted as the amplitude of a singular gauge 
field. Note that this parameter can play an important role for nonstationary problems. A 
similar problem has been studied in [jl2[ but it was wrongly concluded there that the fam- 



ily of such models having eigenfunctions given by Bethe Ansatz coincides with the family 
of models satisfying the Yang-Baxter equation. This point has been already clarified in 

i- 

The family of point interactions for the one dimensional Schrodinger operator — ^ 
can be described by unitary 2x2 matrices via von Neumann formulas for self-adjoint 
extensions of symmetric operators, since the second derivative operator restricted to the 
domain C^(R\ {0}) has deficiency indices (2, 2). In what follows we are going to consider 
only the self-adjoint nonseparated extensions that cannot be presented as an orthogonal 
sum of two self-adjoint operators acting in L2(— oo, 0] and L2[0, oo). The boundary con- 
ditions describing the self-adjoint extensions have the following form 



where 

ad — bc=l, 6,a,b,c,d & M. (2) 

if{x) is the scalar wave function of two particles with spin and relative coordinate 
X. (|T]) also describes two particles with spin s but without any spin coupling between 
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the particles when they meet (i.e. for x = 0), in this case (p represents any one of the 
components of the wave function. The values = 6 = 0, a = (i=l in (|l]) correspond to 
the case of a positive (resp. negative) 5- function potential for c > (resp. c < 0). For 
general a, 6, c and ci, the properties of the corresponding Hamiltonian systems have been 
studied in detail, see e.g. [0, |, |, 0, |2|. 

For a particle with spin s, the wave function has = 2s + 1 components. Therefore 
two particles with contact interactions have a general boundary condition described in 
the center of mass coordinate system by: 



(3) 



where ip and ip' are n^-dimensional column vectors, A,B,C and D are n"^ x matrices. 
The boundary condition (j^) can include not only the usual contact interaction between 
the particles, but also a spin coupling of the two particles if the matrices A, B, C, D are not 
diagonal. These conditions are similar to those appeared in ^ during the investigation 
of finite rank singular perturbations of differential operators. 

The matrices A, C, and D are subject to restrictions due to the required symmetry 
condition of the Schrodinger operator. In fact we should have, for any -u, f G C°°(iR\{0}), 

=< u'iO+),viO+) >(pn - < uiO+),v'{0+) >(^n (4) 



- < u'{0-),v{0-) >(rn + < u{0-),v'{0-) >(sn= 0. 
/^From (H) and (|) we get the following conditions: 

A^D - C^B = 1, B^D = D^B, A^C = C^A, (5) 

where f stands for the conjugate and transpose. Obviously (|1]) is the special case of (|^) 
at s = 0. 

In the following we study quantum systems with contact interactions described by 
the boundary condition (^, in particular, A^-body systems with 5-interactions. We first 
consider two spin-s particles with (5-interactions. The Hamiltonian is then of the form 

H=i-^^-^Jh + 2h6ix,-x,), (6) 



where I2 is the x identity matrix, h is an x Hermitian matrix. If the matrix 
h is proportional to the unit matrix I2, then H is reduced to the usual two-particle 
Hamiltonian with contact interactions but no spin coupling. 

Let Cat a = 1, ...,n, be the basis (column) vector with the a-th component as 1 and 
the rest components 0. The wave function of the system (|^) is of the form 

n 

0a/3(a;i,a;2)e„ (g) 6/3. (7) 

In the center of mass coordinate system, X = {xi+ a;2)/2, x = xi — X2, the operator @ 
has the form 

The functions (j) = (j){x, X) from the domain of this operator satisfy the following boundary 
condition at x = 0, 

n 

^'^^{0+,X)-^'^^{0-,X)= Kx,o.pMO,X), 0„/3(O+,X) = 0,/3(O-,X), a,P = l,...,n, 

a,/3=l 

(9) 

where the indices of the matrix h are arranged as 11, 12, In; 21, 22, 2n; nl, n2, nn. 
@ is a special case of for A = D = I2, B = and C = h. h acts on the basis vector 

n 

of particles 1 and 2 by hca ® = ^ hai3,y\e-y (S> ex. 

7,A=1 

According to the statistics ip is symmetric (resp. antisymmetric) under the interchange 
of the two particles if s is an integer (resp. half integer). Let ki and k2 be the momenta of 
the two particles. In the region xi < X2, in terms of Bethe hypothesis the wave function 
has the following form 

^ = ^i^2e*('=i^i+'^2^2) + U2ie''^'''^'+''''''\ (10) 

where uu and U21 are x 1 column matrices. In the region xi > X2, 

ip = (pi2^i2)e*(*^i^^+^^^i^ + (P^Vi)e^^*^^^^+^i^i\ (11) 

where according to the symmetry or antisjnnmetry conditions, P^^ = p^'^ for bosons and 
P^^ = —p^"^ for fermions, being the operator on the x 1 column that interchanges 
the spins of the two particles. Substituting (|10|) and (0) into the boundary conditions 
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we get 

/ Ui2 + U21 = P^^(mi2 + U21 



\ iki2{u2i - U12) = hP^^{ui2 + U21) + ikuP^'^iuu - M21), ^^'^^ 
where ki2 = {ki — /c2)/2. Eliminating the term P^'^uu from ([12|) we obtain the relation 

U2I = >^21^Ml2 , (13) 

where 

Fa'i' = [2^ku - h]-'[2ikuP^^ + h]. (14) 
For a system of N identical particles with 5-interactions, the Hamiltonian is given by 

TV ^2 N 

^= -E nZ2liV + EM(^*-3;j), (15) 
i=l ^-^i i<j 

where Iat is the x identity matrix, hij is an operator acting on the i-th and j-th bases 
as h and the rest as identity, e.g., hi2 = /i Cg) I3 ® ■■■In, with Ij the n x n identity matrix 
acting on the i-th basis. The wave function in a given region, say xi < X2 < ■■■ < xn, is 
of the form 

n 

^ = </'ai,...,aiv(3;i, •••,a^Ar)eai ® ... e„jv 

ai,...,aiv=l 

= U12 ^e^(''i^i+'=2^2+...+fcivxjv) _^ ^^^^k2Xl+k^x2+...+kt,xt,) _^ (^jyi _ 3) other terms, 

(16) 

where kj, j = 1, ■■■,N, are the momentum of the j-th particle, u are x 1 matrices. 
The wave functions in the other regions are determined from ([T6|) by the requirement 
of symmetry (for bosons) or antisymmetry (for fermions). Along any plane Xi = Xj+i, 
2 G 1, 2, — 1, we have 

where 

y^lt^a, = [2^^a.+i - h,,+,]-'[2jk^^^^^,P^^^' + h,,^,]. (18) 

Here k^-aj+i = {ka^ — ^aj+i)/2 play the role of momenta and P^^+^ = pii+^ for bosons and 
pjj+i — _pjj+i foj^ fermions, where p>^~^^ is the operator on the x 1 column u that 
interchanges the spins of particles j and j + 1. 
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For consistency Y must satisfy the Yang-Baxter equation with spectral parameter 



30, 31 



or 



-i^m,m+l-i^m+l,m+2-iAm,m+l -iAm+l,m+2-iAm,m+lTAm+l,m+2 

ij kj ki ki kj ij i 



-ymr-yrs-ymr -\^rs-\^inr-yrs 

ij kj ki ki kj ij 



(19) 



if m,r,s are all unequal, and 



ij ji ~ ' ij kl — ^kl ^ij 



(20) 



if m, r, s, q are all unequal. By a straightforward calculation it can be shown that the 
operator Y given by (|18D satisfies all the Yang-Baxter relations if 



0. 



(21) 



Therefore if the Hamiltonian operators for the spin coupling commute with the spin 
permutation operator, the A^-body quantum system (0) can be exactly solved. The 

N 

wave function is then given by ([T6|) and ( |17|) with the energy E = y^^kf. 

i=l 

For the case of spin-^, a Hermitian matrix satisfying ( ^1]) is generally of the form 



h'2 



fa ei ei c \ 

el f g 62 

el 9 f 62 

\c* el el b / 



(22) 



where a, b, c, /, ei, 62 G C, (7 G M. We recall that for a complex vector space V, a matrix 
R taking values in Endc{V ^ V) is called a solution of the Yang-Baxter equation without 
spectral parameters, if it satisfies 



'^12'^13'^23 



'7^23'^13'^12, 



(23) 



where TZij denotes the matrix on the complex vector space V ^ V ® V, acting as R on 
the 2-th and the j-th components and as identity on the other components. When ^ is a 
two dimensional complex space, the solutions of (^) include the ones such as Rq which 
gives rise to the quantum algebra SUq{2) and the integrable Heisenberg spin-i chain 
models such as the XXZ model {R corresponds to the spin coupling operator between 
the nearest neighbor spins in Heisenberg spin chain models) |TD|, |T^, ^ . Nevertheless 



in general h2 does not satisfy the Yang-Baxter equation without spectral parameters: 
i_ 1 1 ill 

^12^13^23 7^ ^23^13^12- But (p2D iucludes the Yang-Baxter solutions, such as Rg, that gives 
integrable spin chain models (for an extensive investigation of the Yang-Baxter solutions 
see [13, 16|). Therefore for an A^-body system to be integrable, the spin coupling in 
the contact interaction ([T5|) is allowed to be more general than the spin coupling in a 
Heisenberg spin chain model with nearest neighbors interactions. 

We now investigate the problem of bound states. For N = 2, from (^) the bound 
states have the form, 

^^ = ««e^l---l, « = l,...,n^ (24) 

where Ua is the common a-th eigenvector of h and P^^, with eigenvalue A^, s.t. hua = 
AaUa and c + aA^ < 0, P^^Ua = Ua- The eigenvalue of the Hamitonian H corresponding 
to the bound state (plD is — (c + aK^Y /2. We remark that, whereas for the case of the 
boundary condition (|l]), for a 5 interaction one has a unique bound state, here we have 
11? bound states. 

By generalization we get the bound state for the A^-particle system, 

V'a =^ae"^^'>^'^'"^^', a = l,...,n^ (25) 

where Va is the wave function of the spin part. 

It can be checked that satisfy the boundary condition (^ at Xi = Xj for any 
i j E 1, N. The spin wave function v here satisfies P^-'Va = Va and hijVa = A^fa, for 
any i ^ j. 

It is worth mentioning that ip^ is of the form (|1^) in each of the above regions. For 



instance comparing with ([T6|) in the region Xi < X2... < o^at we get 



^ _^ c + oAq ^^ _ ^ _^ = /;;2 + ic, kN = -h, (26) 

for a = 1, ...,71^. The energy of the bound state is 

K = ~^-^±^N{N' - 1). (27) 

Now we pass to the scattering matrix. For real ki < k2 < ...kN, in each coordinate 
region such as Xi < X2 < ...xn, the following term in (0) describes an outgoing wave 

V'o^t = Mi2...ive'(*^^"^+-+'-"-). (28) 
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An incoming wave with the same exponential as (28) is given by 



in the region xn < xn^i < ••• < Xi. From (0) the scattering matrix S defined by 
'4^ out = Sipin is given by 

5*= [X2lX^i...XN\\[X^2Xi2---XN2\--\XN{N-l)\-, (30) 

where Xij = Y^^ P'K 

The scattering matrix S is unitary and symmetric due to the time reversal invariance 
of the interactions. < s'^Sg... 5^15*13152. ..sat > stands for the S matrix element of the 
process from the state (/ciSi, k2S2, k^s^) to the state {kis[, k2S2, k^s'^^). 

The scattering of clusters (bound states) can be discussed in a similar way as in 
For instance for the scattering of a bound state of two particles (xi < X2) on a bound state 
of three particles (0:3 < < x^), the scattering matrix is 5* = [X32X42X52] [X31X41X51] . 

The integrability of many particles system with contact spin coupling interactions 
governed by separated boundary conditions can also be studied. Instead of (^ we need 
to deal with the case 

(/)'(0+) = G+0(O+), 0'(O_) = G-0(O_), (31) 

where are Hermitian matrices. For = G^ = G, G"^ = G, there is a Bethe Ansatz 
solution to dH) with Y^:^^ in given by 

Let r be the set of eigenvalues of G. For any Aq, G F such that Aq, < 0, there are 
2N{N-i)/2 ]3Q^n(j states for the A^-particle system, 

i'l = Vae W^^k - xi) + tue{xi - x,))e^"^>>. I"'""^!, (33) 

k>l 

where v^e is the spin wave function and e = {eui ■ k > l}; eu = ±, labels the 2^^^"^)/^- 
fold degeneracy. The spin wave function v here satisfies P^^Vae = ^ijVae for any i j, that 
is, p^^Vae = ^ijVat foi bosous and p^^Vae = —^ij^at fof fcrmious. 



Again ip^^ is of the form (|T^) in each of the regions x,^ < < ... < Xj^. For 
instance comparing ip^^ with (|l^) in the region Xi < X2... < Xn we get ki = iXa{N — 1), 
k2 = ki — 2i\a, ks = k2 — 2iXa,...,k]s[ = —ki. The energy of the bound state ipae is 

Ea = -^N{N' - 1). (34) 

We have investigated the integrable models of A^-body systems with contact spin cou- 
phng interactions. Without taking into account the spin couphng, the boundary condition 
(|I]) is characterized by four parameters (separated boundary conditions are a special lim- 
iting case of these). Obviously the general boundary condition (|^) we considered in this 
article has much more parameters. The classification of the dynamic operators asso- 
ciated with different parameter regions is a big challenge. As we have seen, the case 
A = D = I2, B = 0, C = h corresponds to a Hamiltonian with 5-interactions of the 
form (^ (for A^ = 2). It can be further shown that (for N = 2) the following boundary 
condition 

f ip \ f I B \ f ij \ 



^/o+ V /\^/o- 



(35) 



corresponds to a Hamiltonian H of the form: 

H = -Dl{l + B6)-BDJ', 

where B is an ra^ x Hermitian matrix, is defined by (D^/)((/9) = —f(-^Lp), for 
/ G C^{IR/ {0}) and (f a test function with a possible discontinuity at the origin. 

The boundary condition 



describes the Hamiltonian 



H = -Dl + iB{2DJ -S'). 



(36) 



We have introduced boundary conditions depending on the spin states of the particles 
and studied several special cases. A complete investigation of integrable A^-body systems 
and Hamiltonian operators corresponding to the general boundary conditions of the form 
(^ still remains to be done. 
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